Josephson probe of Fulde-Ferrell-Larkin-Ovchinnikov phases in atomic Fermi gases 
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We study theoretically two spatially separate one-dimensional atomic Fermi gases in a double- well 
trap. A Josephson junction between a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superfluid and a 
Bardeen-Cooper-Schrieffer superfluid could be realized by creating a weak link at a given position 
and by tuning independently the spin polarization of gases in each well. Owing to the proportionality 
between the Josephson current and order parameters, the spatial inhomogeneity of the FFLO order 
parameter can be directly revealed through the current measurement at different positions of the 
link. Thus, the Josephson effect provides deflnitive evidence for the existence of FFLO phases. 
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Strongly attractive Fermi gases with imbalanced spin 
components are ubiquitous systems in diverse fields of 
physics [H. They are building-blocks of atomic nuclei, the 
matter in neutron stars and even the quark-gluon plasma 
that comprised the early Universe. Polarized Fermi gases 
also appear in soHd-state superconductors subjected to 
either an internal exchange field or external magnetic 
field. A recent example attracted intense attentions is 
trapped atomic gases of neutral fermions with unequal 
spin populations Owing to the flexibility in the con- 
trol of the constituents and interaction strengths, atomic 
Fermi gases thus provide the most promising place for 
observing many exotic forms of matter. 

The ground state of polarized Fermi gases remains elu- 
sive The mismatched Fermi surfaces in polarized 
environment cannot guarantee the standard Bardeen- 
Cooper-Schrieffer (BCS) mechanism, which requires a 
pairing of two fermions on the same Fermi surface with 
opposite spins. Various exotic forms of pairing have 
been suggested such as Fulde-Ferrell-Larkin- 

Ovchinnikov (FFLO) state with spatially varying order 
parameters deformed Fermi surface [51], interior gap 
or Sarma superfluidity [6|, and phase separation 

Among these, the FFLO state is of particular inter- 
est, since the Cooper pairs may condensate into a fl- 
nite center-of-mass momentum state. The request for 
the FFLO state has lasted for more than four decades 
in many branches of physics. Only recently experimental 
evidences of its existence have been found in the heavy 
fermion superconductor CeCoIns. The observations in- 
clude the speciflc heat, magnetization, and penetration 
depth measurements [a|. However, deflnitive identiflca- 
tions of the FFLO phases are desirable. These should be 
based on phase-sensitive measurements that can directly 
reveal the spatial variations of the phase of the order pa- 
rameter. One possibility is the Josephson effect 

In this Letter we propose an atomic Josephson effect 
to detect the existence of FFLO phases. Theoretical 
studies suggest that such phase is more favorable in the 



low-dimensional systems. Particularly it dominates in 
the one-dimensional (ID) polarized gases (lol . 11 1. We 



therefore consider two spatially separate ID atomic Fermi 
gases with tunable spin polarizations in a tight double- 
well potential on an atom chip, where the lateral motion 
of fermions is frozen, while axial motion is weakly con- 
flned. A weak link at a speciflc position a;o is easily cre- 
ated by superimposing a narrow dipole dimple potential 
to allow tunneling. Fig. 1 presents a schematic view of 
the conflguration, together with the potential and par- 
ticle density proflles. Its detailed realization will be ad- 
dressed later. 

The underlying physics of our proposal is easily under- 
stood using Ginsburg-Landau (GL) theory [l3|. Assum- 
ing that the order parameters or condensate wave func- 
tions in the left and right wells are described by 5'bcs (x) 
and \1/fflo(2;); respectively, the Josephson current by GL 
theory is 
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where J is the characteristic tunneling energy and is de- 
termined by the small overlap of the condensate wave 
functions along the lateral direction. Aa; is the width 
of the narrow tunneling link. The BCS order parame- 
ter in Eq. ^ is essentially spatially independent, while 
the FFLO one is oscillatory in real space with period 
27r?i/gFFLOj where ^fflo is the center-of-mass momen- 
tum. Thus, provided that Ax ^ 27rfi,/(7FFLO, the mea- 
surement of the maximum Josephson current results di- 
rectly '^FFhoix = xq). By displacing axially the har- 
monic traps and consequently changing the position of 
weak Hnk xq, a series of measurements therefore reveal 
the whole spatial inhomogeneity of the FFLO order pa- 
rameter. As shown in Fig. 2, the simple GL picture is 
verifled by much complicated microscopic calculations. 

Our results are obtained by solving mean-fleld 
BogoHubov-de Gennes (BdG) equations for each well, 
while treating the tunneling through the weak link within 
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Figure 1: (Color online) Schematic view of a proposed atomic 
Josephson junction on an atom chip, (a) The double-well po- 
tential landscape. Two needle-like (BCS and FFLO) super- 
fluids are located in the left and right wells, respectively. Re- 
ducing the double- well barrier at a speciflc position a; = a;o by 
superimposing a narrow dipole dimple potential, the Cooper 
pairs in superfluids can tunnel back and forth between wells, 
leading to a Josephson current, (b) Lateral distributions of 
the potential (thick dot-dashed line) and spin up and down 
particle density proflles (thin solid and dashed lines) away 
from the weak link at xo- The inter- well distance in experi- 
ments would be about lOOfim. The potential with superpo- 
sition of a dimple potential is shown by the thick solid line. 
Then, the particle densities and the order parameters can 
have sm all ove rlaps within the weak link. The length scale 
aho = \/h/muj is around 40/im. (c) Axial profile of the po- 
tential at 2: = 0. 



linear response theory. Our calculations are performed 
specifically for atomic Fermi gases. However, as the 
FFLO physics is a fundamental issue that is of impor- 
tance to many research fields, they can have potential 
implications beyond ultracold atoms. 

We assume that the dipole dimple potential is a small 
perturbation and the resulting weak link at xq does not 
distrub considerably the distributions of the order param- 
eter and particle density profile in each well. The atomic 
Josephson junction in Fig. la is then well described by 
a tunnehng Hamiltonian. By integrating out the lateral 
degree of freedoms in a tight-binding approximation, it 
takes three terms: 

H = Ul+Ht + Ub., (2) 
Ht = Vo [^t. (a^o) ^Ra (xo) + H.c] , (3) 

where u the spin index. The terms TLl and Hr 

are respectively the Hamiltonians for fermions on the left 
and right sides of the junction, and can be expressed 
in terms of operators ipL.R-.a (x). They contains all the 
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Figure 2: (Color online) Underlying physics of the Joseph- 
son effect, as a probe of FFLO phases. With a BCS and a 
FFLO superfluid placed on the left and right sides of the junc- 
tion, respectively, the maximum Josephson current is roughly 
proportional to the FFLO order parameter. The number of 
fermions on each well is A''^ = A''^ — 128, while interaction 
strengths 7l = 7i? = 1.6. The order parameter is normalized 
by the Fermi energy Ep = NLhuj/2, and is shifted upwards 
by an amount of 0.50Ep for clarity. The spin polarization in 
the FFLO superfluid is 0.25. Uivjv is the conductance of a 
corresponding normal junction. 



many-body interactions. Assuming that the width Aa; is 
the smallest length scale, in Ht we approximate all oper- 
ators ip (x) ~ 'ip {xq) and introduce a transfer parameter 
Vb — J Ax. This is valid as far as Ax <C 2'Kh/ q^-p-Lo- 
We shall take a small constant transfer parameter, which 
corresponds to the small overlap of two order parameters. 
The overlap could depend weakly on the position of the 
weak Hnk. However, the assumption of a fixed transfer 
parameter is sufficient to capture the qualitative feature 
of the Josephson effect. 

In each well the ground state of an attractive gas of 
N = -\- Ni fermions with polarization P = {N^ — 
Ni)/N is conveniently determined by using the BdG for- 
malism [l3| that describes the quasiparticle wave func- 
tions Uri{x) and Vri{x) with a contact interaction g (the 
well index is suppressed for clarity), 
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is 



where Ti!^ ^ 

the single particle Hamiltonian under axial harmonic 
trap and Hartree potential. The chemical potentials are 
shifted as n-^^i = /i±(5/i to account for the unequal popu- 
lation N^^i. The order parameter A(a;) and are cal- 
culated by self-consistency equations for the gap, A{x) = 
gid'}Zj^UT^{x)v*{x)f{Er^), and for the densities: n]{x) = 

ll^\un{x)t f{Er^) and 711(2;) = Y.^K{x)\^ fi'E^), 
with f{x) = l/(exp[x/A:Br] + 1) being the Fermi func- 
tion. These must be constrained so that / dxn^ ^ {x) = 
iV| |. We note that the unequal chemical potentials break 
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Figure 3: (Color online) Zero temperature mean field density 
profiles (solid lines) of a polarized gas at A'' = 128 and 7 — 
1.6, compared with the results from Gaudin solutions (dot- 
dashed lines) [lol |. In (b) the oscillations in the spin- up density 
outside the FFLO core are due to Friedel oscillations. 



the time-reversal symmetry. Thus, the sum over energy 
levels is done for all the eigenstates with both positive 
and negative energies i?^. Generically, the interaction 
strength is parameterized by a dimensionless coupling 
constant 7 — —mgid/{fi'^no), where ng is the center den- 
sity of an ideal gas. In the weak or intermediate coupHng 
regimes {i.e., 7 < 10°), the mean-field BdG theory ap- 
pears to be very accurate. Fig. 3 shows the mean-field 
results of density profiles of a gas at 7 = 1.6 and P — 0, 
0.25, compared with that from exact Gaudin solutions 
and local density approximation The agreement is 
excellent. In case of a nonzero spin polarization, both 
theories predict a two-shell structure with a partially po- 
larized superfiuid at the trap center and a fully polarized 
normal state at edge. The mean-field order parameters 
are given in Fig. 2. The spatial variation identifies clearly 
that the partially polarized phase at center is indeed a 
FFLO superfiuid. 

The main observable of interest, the rate of trans- 
ferred atoms from, e.g., right well to left well, is de- 
fined by /(<) =< dNL{t)/dt >. In analogy to super- 
conductors where the fiow of electrons out of the super- 
conductor establishes an electrical current, we call / the 
current. By rewriting the transfer Hamiltonian Ht = 
J2cri'^'^ + ^t) where = Tp+^ (xq) ipRo- (xo), the equa- 
tion of motion leads to dNL{t)/dt = i[HT{t), NLit)] = 
i^^[Ao-(i) — A+(t)]. Bearing in mind that the link 
at xq is weak so that the transfer of atoms can be 
treated as a perturbation, we use the linear response 
theory in which the current I{t) — —i dt' < 

[dNL(t)ldt,nT{t')] >0= Haa'^-oo^t' < W " 

A+{t),A^'{f) + A+,{t')] >o. Here the subscript in 
the average refers to the unperturbed systems Hl and 
Hr. Two contributions can be easily identified: the 
normal single-particle current and the Josephson cur- 
rent of Cooper pairs. The explicit expression of the lat- 



ter is given by pj| Ij{t) = exp[i(<^fl - <^l) -|- i2(^fl - 
^^L)t/h] dti ex.Y>[-i{^^R-<y - f^L-a){t - ti)/h] < 

[Aa-{t), A-cr{ti)] >o +C.C., where we introduce an interac- 
tion representation with respect to TYl and Hr, and rep- 
resent it by a tilde in operators. The global phases of or- 
der parameters, ipL and ipR, are made explicitly in I.j{t). 
Their difference, together with the factor 2{^,r — fiL)t/fi, 
drives the dc and/or ac Josephson currents even at the 
zero chemical potential difference between wells. 

With the help of the Wick theorem, in the statistical 
average of I,j{t) we split the four fermionic field oper- 
ators. The integration of the average over time ti can 
then be expressed in terms of the retarded correlation 
functions |13l |. 



X 



[f [E,) - f{E,)] 
^'r ~n + E^ -Ej 



(6) 



where the indices i and j refer to, respectively, the energy 
levels in the left and right wells, and the subscripts L 
and R are the indices for wells. We have abbreviated 
u — u{xo) and v — v{xo). Consequently, the Josephson 
current can be calculated as Il3ll. 
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where the maximum current 
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The Josephson current thus oscillates in phase with a 
peak value I™^^. It is clear from the expressions ^ and 
([6]) that correlation functions become roughly a product 
of two order parameters if we approximate the gap equa- 
tion A(x) cx '^^Uri{x)v'^(x). This particular structure 
emphasizes the microscopic origin of the GL theory of 
Josephson effect. 

A simple expression for J™^'' can be derived when we 
assume identical uniform BCS superfiuids in both wells, 
= (7r/2)awwA, where ^nn = iV;^ ln)[2m/ {fV^Ep)] 
may be viewed as the conductance of a normal junction, 
with Ep being the Fermi energy. Therefore the value 
of gap can be determined by measuring I™^^. In the 
presence of traps or inhomogeneous superfiuids, one has 
to resort to the numerical calculations. Fig. 2 provides 
the most promising scheme, where a BCS superfiuid is 
set in the left well as a reference system, while the order 
parameter of a FFLO superfiuid in the right well is to 
be determined. As a result of the fiat distribution of 
BCS order parameter, we anticipate J™'^'^ oc Afflo(2;o); 
which is confirmed numerically. Therefore, the spatial 
inhomogeneity of FFLO phases can be precisely detected 
by Josephson effect at varying positions of the weak Hnk. 

We are now in position to discuss the experimental re- 
alization of the atomic Josephson effect. The tunneling 
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junction of two ID Fermi gases sketched in Fig. la is 
most easily set up on an atom chip using a combination 
of static and radio-frequency magnetic fields that forms 
an adiabatic double well potential. This technique has al- 
ready been applied successfully to create a matter-wave 
interferometry of ID Bose gases A narrow weak 

link between wells might be built up by the addition of a 
tight dipole dimple microtrap at a specific position. 
By tilting the bias field to displace the center of the axial 
traps, the position of the link can be effectively moved. 
Further, the inter-well distance ~ 100 /im is sufficient 
large to allow an independent adjustment of the spin po- 
larization in each well with a radio-frequency sweep [2|. 

To observe the Josephson oscillations, it is necessary to 
fulfill two conditions: (i) the number of fermions involved 
in the oscillations, to be measured by phase-contrast 
imaging should be larger than the finite optical res- 
olution of the imaging system, but small enough to en- 
sure the validity of the Hnear response theory; and (ii) 
the width of the weak link Aa; should be much smaller 
than the period of FFLO order parameter 27r?i/gFFLO- 
Typically, the total number of atoms in one well would 
be around TV ~ lO'^, and the frequency of axial trap 

^ 27r X 1 Hz, which is much smaller than the radial 
frequency uj± ~ 27r x 10^ Hz so that the ID condition 
Nuj < uj± holds For ^Li atoms, with these param- 
eters we find 27rfi./gFFLO ~ 2TrhvF/^o ^ 30 /xm at an 
interaction coupling 7 = 1.6, where vp and Aq are re- 
spectively the Fermi velocity and BCS gap at the trap 
center. This period is much larger than the width of 
dimple potential Aa; that is about several micrometer, 
and therefore the condition Ax <C 2Trh/qppi^o is well 
satisfied. On the other hand, with a tunneling energy 
J ~ 10~^huj±, the maximum Josephson current is es- 
timated to be /j,max ~ fTATArAo ~ lO'^ s~^. Assuming a 
time scale of oscillation ~ 0.1 s, the number of transferred 
fermions is then around hundreds, an order of magnitude 
smaller than the total number of atoms. This number is 
within the present experimental detection limits, i.e., the 
Josephson oscillation in a Bose-Einstein condensate with 
similar number of atoms and time scale has already been 
observed [lB |. 

Before concluding two remarks are in order concern- 
ing the temperature and interactions, (i) In contrast to 
the 3D case, the ID FFLO state is notably stable in 
response to a nonzero temperature. In magnitude the 
critical temperature of the FFLO state is at the same 
order {i.e., a half or one third) of its unpolarized coun- 
terpart, Tbcs 4.54e~'''/27Tp^ where Tp is the Fermi 
temperature. Given an intermediate interaction 7 — 1.6, 
we estimate a critical temperature Tfflo ~ O.IOTf at 
the trap center, which is well above the lowest tempera- 
ture reported so far [l3| • We anticipate a lower transition 
temperature inside the weak link because of reduced den- 
sity. However, it can be much enhanced by increasing 7. 
(ii) In practise, gid is parameterized by a 3D scatter- 



ing len gth a3rf, gid = 2h'^uj±a3d/ii - Aa3d/a±), where 
a± = y^h/muj^ and A ~ 1.0326. The denominator indi- 
cates a confinement-induced Feshbach resonance that oc- 
curs when a^d o-± as observed experimentally [l3|- 
For ^Li, the background a^d — 114 nm, and a± ~ 130 nm 
corresponding to uj± ^ 2t: x 10^ Hz. Therefore, by finely 
tuning the shape of double well potenital and hence , 
the coupling constant 7 may be varied at will. 

In summary, we have proposed a scheme to realize 
the atomic Josephson junction to detect the existence 
of FFLO phases. Our proposal opens the possibility for 
creating ultracold Fermi gases for practical applications, 
such as precision measurement and interferometry. 
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